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In a recent experiment by Eichmann et ai, polarization-sensitive measurements of the fluores- 

. cence from two four-level ions driven by a linearly polarized laser were made. Depending on the 

^\ • polarization chosen, different degrees of interference were observed. We carry out a theoretical and 

0^ i numerical study of this system, showing that the results can largely be understood by treating the 

, atoms as independent radiators which are synchronized by the phase of the incident laser field. The 

^ ' interference and its loss may be described in terms of the difference between coherent and incoherent 

Q , driving of the various atomic transitions in the steady-state. In the numerical simulations, which 

■ are carried out using the Monte Carlo wave function method, we remove the assumption that the 

' atoms radiate independently and consider the photodetection process in detail. This allows us to see 

00 , the total interference pattern build up from individual photodetections and also to see the effects of 

(N . superfluorescence, which become important when the atomic separation is comparable to an optical 

' wavelength. The results of the calculations are compared with the experiment. We also carry out 

. simulations in the non steady-state regime and discuss the relationship between the visibility of the 

^ 1' interference pattern and which-path considerations. 

g ! PACS numbers; 32.50.-f d, 03.65.Bz, 42.50.-p 

o : 

; 

^ • I. INTRODUCTION 

CJ ■ 

• 1— I . 

Young's two slit experiment is the canonical experiment demonstrating the wave nature of light. With the addition 
(-H , of a measuring device to determine which slit the photon has passed through, the phenomenon of wave-particle duality 
O ' can be investigated. Many such path detection schemes have been proposed in the past, such as the Einstein recoiling 
slit j^] and Feynman's light microscope |Q. 
^ • Recently, interest has arisen in the interference of the fluorescent light from two driven atoms which play the role 
^ ' of the slits in Young's experiment. Various situations have been theoretically studied, using two-level and four-level 
QQ . atoms, and the possibility of placing the atoms in a cavity has also been considered. ||^-|^ 

' In a recent experiment performed by Eichmann et al. the light scattered by two ^^^Hg'^ ions in a trap was 
,, observed. In this experiment, a linearly-polarized traveling wave coherent field was used to drive the ions and the 
' scattered light was observed in two orthogonal polarizations. An interference pattern is observed in one case but not 
in the other. These results were explained in two ways, firstly in terms of a which-path argument based on ways in 
which a single photon can interact with the two atoms and secondly in terms of a theoretical analysis by Polder and 
Schuurmans [Q who considered coherent driving of a single four-level atom. 
^ . In this paper, we carry out an analysis which considers the presence of both atoms and stress that for the coherent 

• ^ ■ steady-state excitation as used in the experiment, the which-path argument is not applicable. In sections || to ^ we 

extend the analysis of Polder and Schuurmans to the experimental situation. In brief, the atomic dipoles are only 
(-H ■ coherently driven in a specific direction given by the polarization of the incident light field, and so interference is only 
expected for scattered light of that polarization. Scattered light of the orthogonal polarization is also present, since 
the atoms can decay spontaneously, but this is incoherent with the driving field and does not lead to interference. 
Besides the above analysis which is carried out using th e ma ster equation, a Monte Carlo wave function simulation 



0\ 



X 



was also carried out, as is discussed in sections VI and VII. We use the method to simulate the system without 



5h ' resorting to the assumptions of large atomic separation and independent radiators used in the previous analytical 



sections. When a fluorescent photon is detected as traveling in a given direction, it may not be possible in general 
to determine from which atom it arose, and this leads to additional "dipole-dipole" type of terms which become 
important when the atomic separation is comparable to the optical wavelength. At small separations, the atoms 
act cooperatively giving rise to superfluorescence whereas for large separations, we smoothly recover the results for 
independent radiators. 

In order to compare our results with the experiment, we include the effects of the classical motion of the ions within 



the trap (Sec. VIII). As expected, this reduces the visibility of the interference patterns, bringing the results more 



closely to those observed. 



In the last section (Sec. IX) we consider a situation where the two atom system is re-prepared to its initial state 



after each individual photon has been detected. This allows us to study the transient regime of this system and we 
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find that interference is present in both polarizations, although one of these disappears in the steady-state. Depending 
on the initial state however, the visibility of this transient pattern varies, and a which-path argument may be used in 
this case to explain the results. 



II. THE MODEL 

The atoms are modeled as four-level atoms (see fig. |l|) interacting with a linearly polarized light field polarized 
along the z axis and traveling along the y direction. We treat the incident laser light as a classical field and the 
external motion of the atoms will be neglected. The separation between the atoms will be considered to be large 
(many wavelengths of the laser light) so that dipole-dipole interactions can be ignored. In this approximation the 
atoms can be considered as independent fixed sources except that their radiation is synchronized by the definite phase 
of the incident laser light. The atomic dipole operator /2 is the sum of atomic raising /x^ and lowering fji^ operators 
whose components are given by 

Mi^/i(|l)(4| + |2)(3|)x (1) 
/.,i = -^M|l>(4|-|2)(3|)y (2) 
Mi = M|2>(4|-|l)(3|)z, (3) 

where is the k component of the atomic dipole, /i is the dipole matrix element, x, y and z are the usual Cartesian 
unit vectors. The atomic operator couples the i'th and j'th levels of the atom. The system Hamiltonian for one 
of the atoms in the interaction picture is 

H,, = ^^lA ((T44 + (T33 - (T22 - CTii) -I- Ti^ ((T24 + 0-42 - - 0-31) , (4) 

where A is the detuning between the atomic transition and the incident laser light. For simplicity, the ground 
states (and the excited states) are considered to be degenerate in energy (see fig. |l|). We define the Rabi frequency 
n = iiEo/21i where Eg is the electric field amplitude of the incident laser and we use the notation Uij = \i){j\. The 
equation of motion for the reduced density operator ps for the single atom is 



dps I 



dt 

where £rciax is the relaxation superoperator 



[Ps,Hs] +£rolax(/5s), (5) 



-Crclax = —7 (f44Ps + PsCTii + CT^sPs + PsC^Ss) 

+7 (cri4/9sCr4i + ((724 - CTls) Ps (0-42 - 0-31) -I- a23PsCr32) , 

and we have assumed that both excited states have the same decay rate 7. 

III. THE OPTICAL BLOCH EQUATIONS 

There are sixteen optical Bloch equations corresponding to the sixteen atomic operators. However, due to the 
geometry chosen for the incident light, the sixteen equations de-couple into two sets of eight. The first set of coupled 
equations involve the populations and the "linear coherences" (((713) and ((724) terms) which relate to the z component 
of the atomic dipole (see eqn. §). 

^A^Tis) = (-7 - *A)(ai3) - li^iM - {<Jii)) (6) 
at 

ii^^si) = (-7 + «A)(a3i) + in{{a33) - (an)) (7) 
at 

3-((T24) - (-7 - «A)(a24) - ini{a22) - {(Jn)) (8) 
at 

^((T42) = (-7 + «A)(a42) + ini{a22) - (0-44)) (9) 
^(cTii) = -ifliiasi) - (ai3)) + 7((a33) + (^744)) (10) 
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^(f722> = -l^iM - (^T42>) + 7(('^33) + (^T44>) (11) 

^(^733) = -*r!((ai3) - (^3i» - 27(a33) (12) 
at 

— ((T44) = -«rj(((T42) - {(J24)) - 27(ct44). (13) 

Since normalization requires that 

4 

1=1 

the number of independent equations may be reduced to seven. Setting the time derivatives to zero for the steady-state 
yields[| 

(15) 
(16) 
(17) 
(18) 

(19) 

(20) 

(21) 

(22) 

where A = A/il and 7 — 7/^1. As expected from the symmetry of the situation, the populations of the two ground 
states are equal as are the populations of the two excited states. 

A graph of the populations versus laser light strength is shown in fig. ||(a). For strong fields the atom is saturated 
with all the populations equal to one quarter. For weak fields, the excited states are close to zero with the ground 
states close to one half. The linear coherences are shown in fig. ^ (b). Since the coherences are purely imaginary we 
plot the imaginary part of (cri3) and ((724) ■ We see that as the field strength increases the coherences grow and reach 
a peak for $1/7^ 2/3. There is a gradual loss of coherence as the field is further increased due to the saturation of the 
atom. At low field strengths the atom spends most of its time in the ground state and is very rarely pumped to the 
excited states. It behaves like a coherently driven linear oscillator and scatters coherently. In the strong field case, the 
atom undergoes Rabi oscillations and is frequently pumped to the excited states from which it spontaneously decays. 
These spontaneous emission events destroy the coherences. 

The second set of equations consist of the "cross coherences" ((cri4) and ((T23)) and the coherences between the 
ground (CT12) and between the excited states (cr34). 

^((714) = (-7 - «A)(ai4> - «0((ai2) + (asi)) (23) 

^(^23) = (-7 - «A)(ct23) + l^{{(T2l) + (fT43» (24) 

^(^12) - + (^32)) - 7(^x34) (25) 

^((734) = -ini{a32) + (fTu)) - 27(a34) (26) 

The other four equations for (1741), (1732), (o'2i) and (CT43) are the Hermitian conjugates of these equations. Solving this 
set of equations in the steady-state we find that only the trivial solution remains so that all these coherences ultimately 
vanish. Any initial coherences involving these terms are damped in the transient regime. Since the cross coherences, 
(cri4) and {(J23), determine the x and y components of the atomic dipole (see eqn. |l|-^, there is no interference from 
light polarized in the xy plane. 



^ These solutions are compatible to those of ref M since only one atom has been included at this time. 
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IV. THE INTERFERENCE PATTERN 



We have calculated in the previous section the steady-state solutions for the atomic coherences and populations for 
the case of a single four-lcvcl atom interacting with a classical laser light field linearly polarized along the z axis. In 
this section we wish to use these to calculate the far-field interference pattern from two such atoms in the case when 
their separation is large enough that they may be treated independently. 

Let us consider a screen placed in the far-field (large y) and oriented in the xz plane. In the far-field, and in the 
paraxial approximation, the electric field at a point on the screen at time t due to an oscillating dipole is proportional 
to the projection onto the plane of the screen of the acceleration of the dipole moment at the retarded time f — r, where 
T is the light travel time from the dipole to the observation point. We shall label a point on the screen by (ti,t2) 
where these are the light travel times from each of the atoms to the point on the screen. The intensity of the light at 
this point is 

I{TuT^)^{ElEi+ElE^:), (27) 

where 

El{t;Ti,T2) cx c-'^^'-^'^^il+c-'^(-*-^-'>ul (28) 

for k G {x, z} , n and U are the atomic dipoles of the first and second atoms respectively and u) is the angular frequency 
of the laser light. If we use the notation o-y for the atomic operators of the first atom and T,ij for those the second, 
we note that the cross terms like {filU^) are zero in the steady-state since 

{nlU^)ss OC ((0-41 + 0-32) (Si4 -I- S23))ss 

= ((0'4l)ss + (0'32)ss)((Si4)ss -|- {T,2s) ss) 

= 0. (29) 

We have factorized the products above since we have assumed that the atoms are independent. The intensity of the 
interference pattern when all the light is detected is given by 

/unpol(Tl,r2) (X +uM + 

+ (fJ-lU^) exp(i[w(ri - T2) - ^o]) 

+ (Ulfxi) exp(-^[a;(rl - rs) - 0o]), (30) 

where the additional phase term (po has been included to allow for the phase of the incident light to differ at the 
positions of the two ions. If this is normalized by twice the intensity due to single atom fluorescence, we obtain 

-fN,unpoi(n, T2) = 1 + Vunpoi cos[a;(ri - T2) - (j)o\, (31) 

where 

K.„poi = I ( j" ] (32) 

We note that this interference pattern is the sum of an incoherent term and a coherent term and so the visibility is 
always less than one half. If we introduce a polarizer at angle i] to the z axis in front of the screen, it is possible to 
separate out the two components. The expression for the intensity is then 

^poi(ti,T2;??) cx (nlfii + UlU^) sin^ T] + (nlni +UlU^) cos^ rj 
+ (ulU^) exp(i[a)(Ti - T2) - 0o]) cos^ t] 

+ {Ulnl) exp(-z[w(Ti - T2) - (ko]) cos^ ri (33) 

Normalizing this as before yields 

-fN,poi(n, T2; ??) = 1 + V^oiin) cos[a;(ri - T2) - (j}o\ (34) 
where the visibility is now given by 
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Thus if we detect the x polarized hght by setting rj = 90°, we see that the visibihty is zero for any field strength 
and detuning. On the other hand if the polarizer is set along the z axis {rj = 0) then the visibility is non-zero. This 
visibility is plotted against field strength in fig. |^ for various detunings with the polarizer aligned with the z axis 
(77 = 0). The maximum theoretical visibility of one occurs for all detunings in the limit of zero field strength. As the 
field strength is increased the visibility falls towards zero, with the detuning determining the rate of this drop. 



V. COMPARISON WITH EXPERIMENTAL PARAMETERS 

In the experiment of Eichmann et al. ^ the incident laser is linearly polarized at an angle of 62° degrees to the z axis 
in the yz plane. This gives a y component to the atomic dipole which is not considered in the above calculation, but 
this does not affect the field on a screen in the xz plane. The inclination also introduces a non-zero phase difference (j)o 
between the two atoms which has be included in our calculations. For more general orientations of the incident field, 
we resort to numerical simulations as described later |^. Another difference between our analysis and the experiment 
is that instead of rotating a polarizer in the output field, a fixed Brewster plate is used and the polarization of the 
incident light is changed. Our analysis is not valid for incident light polarized along the x direction, however our use 
of a rotating polarizer gives some insight into the polarization- sensitive measurement. 

Taking into account the geometry of the experiment, we expect that the measured interference pattern should be 



-^N,unpol ^ ^ ^ 2 



^2 



A +7 



7^2 



+ 7+2, 



cos 



/27rd 



' - 28°) - 



(36) 



when all the light is collected and 



'N,pol 



1 + COS 1] ^ 



,A -^7^ + 2 



/27r(i 



28°) 



(37) 



when a polarization-sensitive measurement is made. Where d is the distance between the atoms, (f) is angle relative 
to the incident light and A is the wavelength of this light. 

The visibility measured in the experiment when all the light was collected had a value of 0.2. The maximum 
visibility from eqn. |3^ is one half. This less than "ideal" measured visibility is not surprising since many detrimental 
factors have not been included in the analytical treatment. Factors such as stray light entering the detector and 
quantum jumps to other levels in the Mercury ions so that we may observe only the fluorescence from one ion. |^ The 
motion of the ions in the trap also smear out the interference pattern, this will be discuss in detail in a later section. 



VI. NUMERICAL SIMULATION OF THE SYSTEM 



In this section, we describe the use of a Monte Carlo wave function method for numerical simulation of the system. 
This approach allows greater flexibility in specifying the geometry of the incident light and also allows us to include 
the photodetection process as part of the simulation. By doing this, we no longer simply relate the far field intensity 
of the fiuoresccnt light to the moments of the atomic variables, but explicitly consider how atomic quantum jumps 
introduce photons into electromagnetic field modes propagating in various directions. The advantage of this approach 
is that it allows the atoms to be separated by arbitrary distances and takes into account the dipole-dipole correlations 
between the atoms. This simulation approach is very closely related to the actual experimental situation as it gives a 
time-resolved classical record of photocounts at different angles which gradually build up into an interference pattern. 

For simplicity, we assume idealized unit efficiency photodetectors in the far field covering the entire 47r steradians 
surrounding the two atoms. Each detector is assumed to be polarization-sensitive and resolves detected photons into 



■^A numerical simulation will be introduced later which can simulate the geometry of the incident light in a more flexible 
manner. 



two orthogonal linear polarizations. Working in spherical polar coordinates with polar angle 9 and azimuthal angle (j) 
(see fig. we use the unit vectors eg and as a basis for the polarization of the emitted photon. These unit vectors 
are related to Cartesian coordinates by 

eg — cos 6 cos (/)x + cos sin ipy — sin 9x (38) 
£0 = — sin 0x + cos (j)y (39) 

where x, y and z are the Cartesian unit vectors. The system Hamiltonian for the two atoms is just the sum of the 
individual Hamiltonians for each atom 

Hsys = ('''44 + ^33 - <^22 - O"!! + S44 + S33 — T,22 — ^n) 

+ TlVlz ((724 + (742 - 0-13 - (731 + S24 + S42 - Si3 - S31) 

+ ((741 + (723 + S41 + 1:23) + Tt-^- ((714 + (732 + ^14 + E32) , (40) 

where we have set — ^^y, and where fl^, Qy and a-re the Rabi frequencies corresponding to the x, y, and 
z components of the electric field 

f^/c = k = x,y,z, (41) 

where fj, is the dipole matrix element and each atom is considered to feel the same electric field E. We are assuming here 
that the incident light is propagating at right angles to the line joining the two atoms. The relaxation superoperator 
is in Lindblad form 

CrelaM = E (^1^™^^ + P^CM + ^ C,nPsCl, (42) 

m m 

where the collapse operators Cm correspond to the couplings between system and the baths. When the atoms 
are far apart, they may be regarded as coupling to independent baths and so collapses associated with one atom 
may be distinguished from those associated with the other. If we explicitly include the photodetection scheme, 
the appropriate collapse operators are parameterized by the direction and polarization of the outgoing photon. For 
spontaneous emission from a single two-level atom, the appropriate collapse operators (including recoil) are |^ 

g \ 1/2 
37^ 



Cn,e = i^-^j exp (-zk • r) (e* • S^) , (43) 

where R is the atomic position operator, k is the wave number of the outgoing photon traveling in direction is 
an operator proportional to the atomic lowering operator (S^ cx <j^) and e labels the outgoing polarization. For our 
situation with two fixed atoms, R may be replaced by the locations of the two atoms fi and fa (see fig. ^) and both 
atoms couple to each mode so that 

Co,, = M {exp (-zk • fi) (e* • Sj) + exp (-zk • fa) (e* • S^) } (44) 

where A/" is a normalization factor. S{ and are proportional to the atomic dipole operators for the first and second 
atoms respectively. The Cartesian components of these operators are 

S{ (X a^x + ayY + cr^z (45) 

cc E,x + S]j,y + E,z (46) 

where the ai's are given by 

(^x (714 + (723 (47) 

(7y = -Z((7l4 - (723) (48) 

(7^ = (724 - (713, (49) 

and the S^'s related to the Sy's in a similar way. We see explicitly from these expressions that when a quantum jump 
(i.e., a photodetection) occurs, this may come from either atom, the coefficient being dependent on the polarization 
and the appropriate phase factor. 
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Between the times of the quantum jumps, the atomic wave function is evolved using the non-Hermitian effective 
Hamihonian given by f\ 



(50) 



On the right-hand side, the summation is responsible for the decrease in the norm of the Monte-Carlo wave function. 
Its expectation value is proportional to the intensity of the fluorescent light and the probability that a quantum jump 
occurs. For our photodetection model, this eventually evaluates to 



/ E 4,.-Co.e~dJ7 = yAA2 {(at a, + at a, + aja. 

e±k 



Sis.) 



-A{a) (at Ea; + crl^y + Sj, a^r + stay) 
^Bia) (atS, + Sta,)} 



where 



A[a) ~ g [~ sma — T(a^ 
Bia) = ^r(a) 



(51) 

(52) 
(53) 



where a ~ |kl • \v2 — fi| is the 27r times the separation between the two atoms in optical wavelengths. The function 
T is given by 



T(a) = — 5- ( — sin a 



cos a 



(54) 



The graph of the functions A and B are shown in fig. ^. Note that we have used the scaled parameter a — ajl-n 
corresponding to units of wavelengths. At a equal to zero where the two atoms coincide, both functions are equal 
to one. As the separation between the atoms is increased both display damped oscillatory behavior decaying to zero 
in the limit of large separations. B is damped more rapidly and is almost zero for separations of more than three 
wavelengths, whereas A has larger amplitude oscillations which can still be clearly seen even after six wavelengths. 
Notice that B corresponds to the size of the correlations between the z components of the atomic dipoles whereas A 
determines the size of the correlations between the x and y components. In the limit of large separation between the 
two atoms, only the first term of eqn. pll remains and we have 



e_Lk 



= (at a, + alay + ala, + ^t + sts^, + Ets^) 

00 

= A/'287r(a44 + a33 + S44 + S33). (55) 
This may be compared with the form expected for two atoms coupling to independent baths for which 

{(J44 + a33 + E44 + E33) 



We see that in the limit of large separation, the atoms behave independently, and that the normalization condition is 
Af — {-f/iny/^. In the opposite limit of the two atoms coinciding 



^see Maimer et al. g 

''The function T(q) has similar structural form to the dipole-dipole term used by T.G. Rudolph et al. [H and J. Guo et al. 



[0 to modify the decay rate of the atoms. 
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^^^Y + ''y"'' + + + ^^^^^ + ^^^^ 



E,|' + |fTj, + Ej,|' + |(7. + I].n . (56) 



This gives four times the decay rate of a single atom instead of two times for the case of large separation since the 
interactions between the atoms grow as the separation decreases until at zero separation the size of these dipole-dipole 
terms become of the same order as the population terms. This gives rise to superfluorescence where the two atoms 
behave as one. This is clearly illustrated by the factorization of the above expression into the coherent sum of the 
two atoms. 

In order to carry out the Monte Carlo simulation, a wave function is evolved using the Schrodinger equation with the 
non-Hermitian effective Hamiltonian until the square of the norm of the wave function reaches a threshold drawn from 
a uniform distribution lying between zero and one. At this time, it is necessary to select one of the possible quantum 
jumps. The probability density of detecting a fluorescence photon traveling in direction (0, 0) with polarization e; is 
given by 



- e OY(j). (57) 



Conventionally, a Monte Carlo simulation requires many runs to be carried out to obtain the desired statistics of 
the system from an ensemble of quantum trajectories. For the problem being considered, the result of one quantum 
trajectory gives the photo-emission record for the two atom system, and we can let this single trajectory run as long 
as it is necessary to obtain the required number of photons. 

VII. NUMERICAL RESULTS 

The raw output of these simulations is a list of the emission times, directions and polarizations of each emitted 
photon. The direction is represented by the two angles 9 and (p i with the polarization aligned either with the unit 
vector €0 or e^. 

The wavelength of the incident laser light used in all of the following numerical simulations was IQAnm. In the 
first simulation we have chosen parameters which gave maximum visibility of the interference pattern. The laser field 
was weak, the ratio of Rabi frequency to decay rate being 0.2. The atomic separation was set at 0.5/zm so that five 
intensity peaks would be observed as 9 is varied through 180°. This gives a sufhcient number of points in order to 
clearly distinguish each peak with the number of photodetections simulated. The incident light is propagating along 
the y axis. 

The intensity pattern as seen on a sphere around the ions is shown in fig. ^ for incident light polarized in the z 
direction. Since the system is symmetric about the z axis, fig. ^ (a) shows an interference pattern as 9 (but not 0) is 
varied. Fig. ^(b) is another simulation in which the incident light is linearly polarized along the x axis. Even though 
the system is no longer symmetric about the z axis, there is little dependence of the interference on the angle cj). We 
shall henceforth only plot the interference patterns as functions of 0, and integrate over the angle (/). This integration 
over (f) modulates the intensity distribution with a sinusoidal envelope since the solid angle integrated over at the 
equator {9 = 90°) is greater that at the poles {9 = 0° or 180°). 

For the case of incident light polarized in the z direction, fig. |^a) shows the intensity distribution for ee polarized 
light as a function of angle 9 from a simulation consisting of one hundred thousand quantum jumps. Out of these, 
about half were detected with eg polarization. The histogram^ is computed using a bin width of 1°. An interference 
pattern is clearly visible. From eqn. |3^, we expect a visibility of 25/27 or about 0.9 but the visibility measured from 
the simulation is approximately 0.7. However eqn. ^is only valid in the paraxial approximation and only considers 
photons polarized in the z direction. When this approximation breaks down (i.e., when 9 is no longer close to 90°), 
the €0 polarization will include contributions from the components of the atomic dipole other than ct^. These will 



'We have plotted the intensity distributions as histograms without the internal lines for clarity. 
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reduce the total visibility as we do not expect interference from the x and y components of the atomic dipole which 
are not coherently driven. If we now look at the orthogonal polarization (the polarization) as shown in see fig. 0(b), 
there is no interference pattern. The polarization is always perpendicular to the z axis for any angle 9 so only the 
X and y components of the atomic dipoles contribute. 

We now consider the case where the incident laser light is polarized along the x axis instead of the z axis. The 
histogram for the intensity of ee polarized light versus 6 is shown in fig. ^a). There is no obvious interference pattern 
now since the 9 polarization consists mainly of the z component of the atomic dipole for angles around 9 — 90°. 
Small peaks can be seen at about 40 and 140 degrees, which are due to the x component of the atomic dipole which 
contribute in larger amounts as the 9 deviates significantly from 90 degrees. Fig. ||(b) displays interference as expected 
since the (j) polarization only consists of the x and y components of the atomic dipole. Since the incident light is 
linearly polarized along the x axis, the y component of the scattered light is not coherent and the visibility should be 
less than seen in ^ (a). From the central peaks the visibility is about one third, approximately half that for the case 
of light polarized in the z direction. 



VIII. EFFECTS OF MOTION 



We have so far ignored the effects of the motion of the ions as this has been secondary to our goal of demonstrating 
how the visibility of the interference pattern varies with polarization. However it is a significant effect in the experiment 
performed by Eichmann et al. We shall model the motion of the ions as harmonic oscillators for each of the three 
main modes, one involving stretching and two involving tilting |^,^. The stretching mode involves motion of the 
atoms along the z axis whereas the two tilting modes involve motion in the xy plane. Thus the stretching mode alters 
the distance between the ions and the period of the interference pattern while the tilting modes change the angle 9 of 
the photon emission displacing the interference pattern. Both of these effects smear the interference pattern, reducing 
the visibility. 

The temperature of the ions in the trap were around the milli-Kelvin region and the frequency of oscillation of 
the modes were of order of a MHz. By comparing the thermal kinetic energy of an ion with the separation 

of the trap energy levels, huj, we find that kBT/{2huj) « 10^ and so the motion may be treated classically. The ions 
experience micromotion within the trap and the recoil momentum of the photon emissions may be neglected. 

In the actual experiment, the values used for the separation between the ions was around 5/im. These values are 
about ten times larger than the values used in the simulations. Thus we expect to see over fifty peaks over a range of 
one hundred and eighty degrees. We carried out a simulation with an ion separation of 3.35/im, where one hundred 
thousand jumps was performed. In fig. |^(a), we show the results for incident light which is linearly polarized along 
the z axis, with no motion of the ions. The peaks are poorly resolved and an excessive length of time would be 
required to obtain sharper resolution for such finely separated peaks. Note that in this graph, the detector is taken 
to be polarization insensitive and we count all the photons. Fig. ^a) shows the correct period (nine peaks from 75 to 
105 degrees) with a visibility of about 0.35 which is lower than the predicted value of 0.46 using eqn. 

The geometry of the incident light in the simulation differs from the that of the experiment. The light is propagating 
along the y axis in the simulation whereas the experiment has the incident light at an angle of 62° relative to the z 
axis. However, the experiment does not detect light over the entire 47r steradians, only a small portion is actually 
observed. This observation region is a small disk perpendicular to the y axis. By summing over the (j) angle in the 
simulation and considering only 9 angles close to the y axis we have an approximation of the experimentally observed 
region of interest. 

The micromotion is simulated in two ways. We keep track of the oscillatory motion of the stretch mode starting 
with some arbitrary initial value. Since we are interested in the long time (steady-state) regime the actual initial value 
we choose is unimportant. The previous numerical procedure is re- used but the atomic separation is altered for every 
photon emission according to the displacement of the stretch mode. This is valid because the characteristic time scale 
of the stretch mode Tstrctch (^ 10~^s) is far longer than the atomic half-life Ti/2 (~ 10~*). The steady-state regime is 
satisfied since the duration of a simulation, Tgim, producing one hundred thousand photodetections is approximately 
10~^s which is much longer than the time scale of the stretch mode. The effects of the tilting motion is to shift 
the interference pattern. Thus during the simulation whenever we calculate a photodetection, the direction of the 
photon needs to be shifted according to the current displacement of the tilting modes. The x and y displacements 
are randomly generated from a harmonic oscillator probability distribution. In hindsight we could have generated z 
displacements from a harmonic oscillator distribution as well, however in the simulations we have chosen to keep track 
of the stretch mode starting from some arbitrary value. Both methods are equivalent for the steady-state regime. 
The requirements on the time scales involved is 

Tsim 3> Tstrctch > Tl/2- (58) 
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As expected, when classical motion is included in the simulation, the visibility of the interference is reduced, see 
fig. 11(b). We see an interference pattern with the peaks placed in the same positions as fig. ^(b) but with a reduced 
visibility of about 0.28. The visibility for individual peaks degrades as one deviates from the central peak at ninety 
degrees. This maybe explained by the fact that we have modeled the motion along the z axis as a stretching where 
the center-of-mass remains stationary, thus the central peak is in the same position whatever the period. As one 
deviates to the sides, a small change in period has a increasingly larger effect on the position of the peak. The 
parameters used in this simulation for the motion were: stretch frequency uj stretch /'^t^ = \/3 MHz, tilt frequencies 
iotiitj'i''^ = V^/"^ MHz (for both modes), radial confinement was 30nm and the amplitude of the stretch mode was 
150nm. Experimental results were shown in the paper by Eichmann et al. for this set of parameters. The visibility 
was about 0.2 which is less than the value obtained in the simulation. The greater visibility in the simulation is not 
unexpected since the effects of stray light and the fact that ^^^Hg'*' ions are not ideal four-level atoms would further 
reduce the visibility. 

IX. THE TRANSIENT REGIME 

So far we have investigated the steady-state regime which corresponds to the experimental situation. However the 
transient regime is also interesting because it is then possible to consider the which-path argument which explains the 
disappearance of an interference pattern as being due to information of which atom has undergone an emission being 
somehow recorded in the system. This knowledge can be obtained for some initial atomic states but not others. 

In this section, we consider that the two atoms are initialized to some known state and then are exposed to the 
incident light. When the first fiuorescence photon is detected, the experiment is stopped and the atoms are re- 
initialized to the original state. Over many trials, an interference pattern may develop in the far field, and it is this 
pattern which is of interest. 

Let us first consider the case where the initial joint atomic state is \ip) |1, 1) and the incident light is polarized in 
the z direction so that the 1^3 transition is driven. When one of the atoms has been excited to the state 3, it may 
decay either to the state 1 or to the state 2. If the decay is to state 1, the final state of the atoms is again |1, 1) and 
it is not possible to tell from which atom the light came from. We thus expect that the detection of light polarized in 
the €0 direction which is sensitive to radiation from the 3^1 transition will exhibit interference. On the other hand, 
if the excited atom decays to state 2, the final state of the atoms is either |2, 1) or |1, 2), leaving a record of the atom 
responsible for the emission. We thus expect light detected with polarization in the direction to not exhibit any 
interference at all. 

The results of a simulation with this initial state are shown in fig. p^(a) for ee polarization and in fig. ^(b) for 
£0 polarization. The atomic separation was 500nm and the ratio of Rabi frequency to decay rate is 0.2 as before. 
As expected, a very distinct interference pattern can be seen in (a) consisting of five peaks. The visibility is about 
0.88, this is quite close to the expected theoretical value of about 0.9. For the same set of parameters this transient 
simulation gives higher visibility than the corresponding steady-state simulation. Also as expected, no interference is 
seen in (b). 

Let us now consider the initial joint atomic state — (|1, 1) + |1,2) -I- |2, 1) 4- |2,2)) /2. In this case, it is not 
possible to distinguish the atom responsible for emitting the photon whether it is detected in the ee or the 
polarizations. We thus expect to see interference for both polarizations. The result of the simulation for the 
polarization is shown in fig. |ll|, confirming that this interference does occur in the transient regime. As discussed 
previously however, since the driving field is polarized in the z direction, this coherence must approach zero in the 
steady-state and there can be no steady-state interference pattern, although interference does occur in the transient 
regime. 

We have used extreme parameters so far, for the case of the |1, 1) initial state we know for certain that each atom 
is initially in a particular ground state and for the other initial state (1-0) = (|1, 1) + |1, 2) -I- |2, 1) -I- |2, 2)) /2) we have 
no initial preference for which ground state each of the atoms are in. In fig. |l^ results for the polarization are 
shown for an initial state lip) = (2|1, 1) -I- |1, 2) -I- |2, 1) -I- |2, 2)) In this intermediate case, there is some visibility 
in the transient interference pattern, but this is only about a third, lower than that seen in the previous case. 

In the transient regime we expect interference in the undriven polarization which does not appear in the steady-state 
regime for the majority of the possible initial states. The interference we see here is due to the presence of initial 
coherences. It is similar to the case of two two-level atoms prepared in the excited states, one would eventually decay 
emitting a photon in the process. The atoms are then re- prepared in the excited states and the process is repeated 
to build up an intensity pattern. It is not important in this regime whether there is coherent driving or not but only 
that the atoms are excited to higher levels thus giving them an opportunity to undergo spontaneous emission. 
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X. SUMMARY 



We have analyzed the experiment by Eichmann et al. and presented theoretical models for understanding the 
degree of interference observed for different incident and detected polarizations. For a simple geometry of the incident 
field, it is possible to solve the optical Bloch equations analytically for a single four-level atom in the steady-state, 
and to extend this result to the case of two atoms, assuming that they radiate independently. Prom this analysis, the 
fluorescent light is seen to consist of a coherent component which exhibits interference and an incoherent component 
which does not. These components may be separated out by polarizers before the light is detected. 

In order to model the actual geometry of the source polarization used in the experiment, a numerical technique 
based on the Monte Carlo wave function method was used. In this method, the photodetection process is simulated 
so that the interference pattern builds up gradually from a series of detections at various angles. This qualitatively 
reproduces the experimental results and also allowed us to consider dipole-dipole interactions which become important 
when the atomic separation is reduced. 

To more closely simulate the experimental configuration, classical motion of the ions was included in the simulation. 
This reduces the visibility of the interference to a value closer to the experimental result. The residual discrepancy 
is probably due to the approximation of modeling the ions as four-level atoms and experimental factors such as the 
presence of stray light on the observation screen. 

Finally, we demonstrated the formation of a transient interference pattern in situations where there is no steady- 
state pattern by resetting the atomic state to a specified initial condition after each photodetection. In this case, a 
which-path argument can be used to understand the visibility of the transient patterns. 
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FIG. 1. Schematic diagram of the four- level atom with decay rate 7 and detuning A from incident laser light of Rabi 
frequency f2. 

FIG. 2. (a) The atomic populations versus field strength. The (an) and (0-22) curves arc graphed as the solid line (they 
coincide) with the ((T33) and ((T44) curves graphed as the dashed line (again they coincide), (b) The atomic linear coherences 
versus field strength. The (0-13) and (a42) curves are graphed as the solid line (they coincide) with the (asi) and (0-24) curves 
graphed as the dashed line (again they coincide). Detuning is set at zero for both graphs. 
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FIG. 3. Visibility plotted against field strength for three detunings. We have the axis of the polarizer aligned with the z 
axis so that the angle 77 is equal to zero. The solid line (a) is the case of zero detuning, the dashed line (b) is A = 67 and the 
dashed-dotted line (c) is for A = IO7. 

FIG. 4. The spherical coordinate system used in the Monte Carlo simulations with the te and e<^ unit vectors defined. We 
also have atoms one and two at positions ri and ri respectively. 

FIG. 5. The function A is shown by the solid line. The dashed line is the plot of the function B and the parameter a is the 
number of wavelengths (of the incident light) separating the two atoms. 

FIG. 6. Plots of intensity as a function of both angles Q and 0. The incident light is linearly polarized along the z axis for 
(a) and the x axis for (b). These plots represents the light pattern distributed over the entire sphere around the two atoms. 

FIG. 7. Histogram of the number of detected photons with (a) d and (b) polarizations as a function of the angle Q. The 
incident laser light was linearly polarized along the z axis. The detuning was set to zero, the atomic separation was 0.5 /im 
and the decay rate 7 = 50. The number of photons detected was 50, 086 and 49, 914 for the B and polarizations respectively. 

FIG. 8. Histogram of the number of detected photons with (a) B and (b) 4> polarizations as a function of the angle Q. The 
incident laser light was linearly polarized along the x axis. The detuning wEis set to zero, the atomic separation was 0.5 /xm 
and the decay rate 7 = 50. The number of photons detected was 49, 770 and 50, 230 for the Q and </> polarizations respectively. 

FIG. 9. Histogram of the number of detected photons irrespective of there polarizations as a function of B (a). Classical 
motion is included for (b). Both simulations used atomic separations of 3.35 /um and decay rates of 7 = 5f2. 

FIG. 10. (a) The Q and (b) ^ polarized light distribution versus the angle d. We have calculated 100, 000 quantum trajectories 
in this transient simulation using the initial state = 1) + |1,2) + |2, 1) + |2,2))/2. There were 49,988 and 50,012 6» and 
4> polarized photons detected respectively. 

FIG. 11. The (i) polarized light distribution versus the angle Q from the transient simulation using the initial state = |1, 1). 
There were 50, 335 </> polarized photons from the 100, 000 quantum trajectories. 



FIG. 12. The 4> polarized light distribution versus the angle 6 from the transient simulation using the initial state 
lip) = (2|1, 1) + |1, 2) + |2, 1) + |2, 2)) /v^7. There were 50, 175 polarized photons from the 100, 000 quantum trajectories. 
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